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ANALYTICAL AND NUMERICAL INVESTIGATION OF STRUCTURAL 
RESPONSE OF COMPLIANT WALL MATERIALS 

By 

R. Balasubramanian 1 


SUMMARY 


Theoretical analysis of an electrostatically driven wall system for a 
compliant wall drag reduction program is reported. The electrostatic wall 
system is capable of producing deflections of many orders greater than the wall 
thicknesses and at small wavelengths. An intermediate large response- theory 
is used for structural analysis. The theoretical predictions are compared to 
bench test results, and good agreement between the two is obtained. The effects 
of aerodynamic loads and perturbation electric fields on the theoretical 
solutions are considered. It is shown that for very small wavelengths 
(X « 2 mm) the aerodynamic effects can be estimated using potential theory 
without loss of accuracy, and the perturbation electric fields do not affect 
solutions as long as the deflections are less than one percent of the wave- 
length. Resonance effects for this type of structure are shown to be fairly 
small. 


1 . INTRODUCTION 


Details of the compliant wall drag reduction program at Langley have 
been discussed in a supplementary report (ref. 1) under the present grant 
NSG 1236. It was pointed out in that report that passive walls with short 
wavelengths and large amplitudes are extremely difficult to design. Extension 
of grant NSG 1236 was given in order to design controlled active wall experi- 
ments. The amplitude of surface motion desired was given to be in the range 
of 5 x 10“ 5 m to 2 x lO" 1 * m, the wavelength of the surface motion in the range 
2.5 x 10” 3 m- to 6 x 10" 3 m, and the frequency range to be between 300 Hz and 
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2 kHz. After careful evaluation of existing techniques for active wall 
experiments (refs. 2 to 4) it was decided to develop an active wall system using 
electrostatic forces as loading to the structure. The choice of materials 
for this system was narrowed to extremely thin elastomers, and the amplitude 
constraint given above dictated large values of the ratio of amplitude thick- 
ness. A nonlinear structural response analysis was conducted to determine 
accurate surface motion predictions. In sections 2 to 4 the electrostatic 
wall system is discussed in detail, and in section 5 comparisons with experi- 
mental measurements are reported. 

2. ANALYSIS OF THE ELECTROSTATIC WALL PROBLEM 

The recent summary paper (ref. 5) suggests that low-speed air experiments 
for compliant wall drag reduction be conducted with controlled or active wall sur- 
face motion to assess the nature of possible turbulent boundary layer modifications 
due to the wall motion. Previous active wall experiments reported in the litera- 
ture (refs. 2 to 4) used mechanical drivers; these drive systems are inadequate for 
producing high-frequency, short-wavelength motion contemplated for active wall experi' 
ments at Langley. The present work describes an electrostatic wall designed to 
operate in a frequency range of 200 Hz to 10 kHz with two-dimensional standing waves 
of wavelength 2 x 10“ 3 to 10“ 2 m. The structural surface is basically a thin elec- 
trically conducting elastomer membrane with a series of transverse electrodes etched 
on a PC board as exciters . The structure is periodic , supported at discrete lines 
by transverse ribs. Figure 1 shows the electrostatic wall system along with the 
electrical hookup. Referring to figure 1, the output from the transformer T 
is biased at the center-tap and connected to the terminals A and B as indicated. 

The conducting membrane surface of width b and thickness h and isotropic 
properties (Young's modulus E, density p, Poisson's ratio v) is supported 
structurally at separations of length l, where i « b; the membrane surface 
is electrically grounded. The electrodes are equally spaced from one another 
and at a separation H from the membrane surface, and they are connected 
alternately to terminals A or B. Each periodic bay of the conducting surface 
has a sealed cavity of volume V(H x £ x b) underneath it. 

The electrostatic wall model is designed to operate at a frequency range of 
200 Hz to 10 kHz. The largest dimension of the model is about 0.4 m, speci- 
fically for testing at the 7 in. x n in. tunnel facility at Langley. The 
largest nondimens ional speed is 
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WL 0.4 x 2-rr x 10000 


« 1 


( 2 . 1 ) 


o 3 x 10 8 

where oj = 2irf and c is the speed of light. Hence, the electric field 
max o 

between the electrodes and the membrane is quasi-static. 

When the electrodes carry voltages, an electric field is set up between 
the membrane and electrodes and the membrane is subjected to a force field. 

The membrane deflects under this loading thereby altering electric field 
distribution. For the case where the electric field is only a weak function 
of the surface undulation, it is possible to uncouple the electric field into 
a primary field (field with no structural motion) and a perturbation secondary 
field which is dependent on the amplitude of the structural motion; The deflec- 
tion of the surface can also be split into a primary deflection (under the loading 

due to the primary field) and a perturbation field, i.e., 

E = Eq + Ex (w,Eq) 

w = w (E q ) + Wi(w(Eq), E) (2.2) 

where Ei << Eo and hence wi << to(Eo) 

3. STRUCTURAL ANALYSIS OF THE MEMBRANE ON PERIODIC SUPPORTS 


3.1. Basic Approximations to Structural Analysis 

We shall make the following assumptions with regard to the structure under 
consideration : 

(i) The structure is a thin elastic membrane with isotropic properties 
(Young's modulus E, density p and Poisson's ratio = v) • 

(ii) The structure is rectangular, flat, and simply supported periodically 
at distances l. Each periodic bay is identical with regard to the loading 
on it, etc. 

(iii) Beneath each bay of the structure is a cavity which is filled with 
an incompressible fluid. Hence any transverse motion of the structure should 
be so as not to decrease the overall volume of this cavity. 
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3.2. Justification of the Basic Approximations 

We shall now examine the above assumptions for the case of a typical 
electrostatic wall as designed; 

(a) The elastomer used was a membrane of thickness between 2.5 ym to 25 ym. 
(Four thicknesses were actually used: 25 ym, 12.5 ym, 6.25 ym, and 2.5 ym) . 

The membrane was uniform and was an elastic material with v = 0 , 3 , p = 148 N/m 3 , 
and E = 4 x 10 8 N/m 2 . In order to make the membrane conducting, it was 
aluminized on one side. The process of aluminizing did not change its uniformity 
nor its elastic properties in any considerable way, hence the validity of the 
assumption (i) . 

(b) The membrane was usually stretched smooth and placed flat over supports 
which were nylon threads. The supports were spaced equally from one another. 

The nylon threads were cylindrical in shape and were glued to the membrane with 
a uniform coating of epoxy resin of negligible thickness. In all cases of 
construction the flatness of the surface was checked using an optical setup and 
was found to be extremely good. The excitation field on the membrane was 
obtained using a symmetrically arranged array of electrodes which were etched 
on a PC board. The electrodes had a flat geometry, were of identical thick- 
ness and breadth, and were equally spaced. The alternate electrodes were 
connected to terminals A and B of a voltage source. Negligible current was 
drawn by the resistance of the wires, and hence a constant potential difference 
existed between terminals A and B. The electrostatic loading on the membrane 
was therefore identical between bays for the form of applied voltages and 

V (see section 4 for the analysis of the electrostatic loading). Other 
forms of loading that might occur on the bays include fluid loadings due to 
static pressure differential appearing across the membrane (the static pressure 
difference can be held constant or nulled using control valves for adjusting 
the back pressure of the cavities) and dynamic loadings induced by turbulence 
(if tested in a wind tunnel with flow over the membrane) and loadings due to 
fluid structure interaction. The main body of the present analysis and 
experimental verification was for bench models where the fluid loadings were 
zero, and hence the fluid loading effects are not considered here. However, in 
a later section we include these cases and indicate how such cases can be 
included in analysis for an accurate prediction of the ensuing structural motion. 
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(c) The cavity underneath the membrane was of uniform depth of about 
6 urn or 10 um. For a membrane with a width of about 0.2 m and periodic 
length of 0.01 m or less, the volume of the cavity is 2 x 10 -8 m 3 for the 
worst case. For a volume change of 0.01 percent of the total volume by any 
form of motion of the structure, the bay would have to suffer a loading 
of 1 kg/m . Since the electrostatic forces acting on the membrane in each bay 
are orders of magnitudes lower in strength, such a change in volume is not 
possible during the motion of the membrane. In other words, the spring stiff- 
ness of the cavity for motions which tend to change the volume is many orders 
of magnitude greater than the stiffness (of the structure + cavity) for 
motions where no such compression occurs (i.e. the stiffness of the structure 
in a deep cavity or in vacuo). For the example considered above, the stiffness 
introduced by the cavity for motions where change in volume occurs is 


k x = AF/AV = 10 9 kg/m 3 , 

which is very large compared to the in vacuo stiffness of the structure 
Based on the above discussion it becomes apparent that assumptions (i) , (ii) , 
and (iii) in section 2.1 are justifiable approximations to a structural 
analysis. 


3.3. Structural Response Theory for Transverse Motions 
of the Periodic Membrane 

The undamped structural response in transverse motion of a simply supported 
rectangular periodic membrane is given by 



(3.1) 


where the quantity is the flexural rigidity of the membrane, p e is 

12(1 - v 2 ) 

the total loading (in vacuo external load and the fluid loading) , N , N , N 

x y xy 

are the midplane forces, and w is the transverse motion of the membrane. 

To evaluate the midplane forces, one must have an idea of the nature of the 
loading on the structure. Based on a classification of the external loading we 
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can classify three distinct regions of structural response. To fix the above 
notion , let us use the following nondimensionalization. 



For the case where the initial stretching forces applied at the edges of 

the membrane (i.e. the forces N , N , N are all positive) the effect 

xo yo xyo * ' 

of these tensile forces is to limit the amplitude of motion and thus raise the 

effective stiffness to transverse motion of the structure. When the initial 

fields N , N , N are negative, the fields are compressive and instability 
xo yo xyo 

of the structure may occur (i.e, buckling instability). 

When the initial compressive stresses are large enough to cause buckling, 
the rectangular membrane will no longer be flat, but will have a deformed shape 
under these compressive loads. With the presence of the cavity the initial 
shape of the membrane is very complicated possibly with local buckling within 
the bay at some points, etc. This situation is rather unpredictable because 
of precise knowledge of the initial tension or compression, especially for 
such flimsy structures. Since it is very difficult to apply compressive 
stresses while mounting, usually care is taken to keep the membrane with 
zero initial tension. The subsequent analysis and prediction are simplified 
by the assumption of zero initial tension. 

The following classification of structural response requirements is made 

with the assumption that the periodic bay is short compared to the width and 

the initial tension field is zero, i.e., b >> l and N = N = N =0. 

xo yo xyo 
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Regime I: small amplitude theory, - The region of applicability of this 


theory is for — = a/y <Q a lG; thus, this analysis is valid for load parameter 
h 


P < 
r e 


1.6tt 4 


12(1 - v 2 ) Y' 


— . The structural response can be conducted using N = N = 


N = 0 in equation (3.3). The transverse motion is therefore decoupled from 
xy 

inplane motions for this case. 


Regime II: (intermediate) large response theory. - For this regime the 

midplane forces are no longer negligible. For analysis of this case, 

Von Kantian' s theory is employed. For the type of flimsy materials we use, 
even at very large amplitudes the structure is within elastic limits. Based 

on an order of magnitude analysis of terms such as u 2 compared to terms like 
2 

u , w we suggest that the intermediate large response theory we develop is 
x x 0 03 

valid for ^ — or for load parameter p g = 0.008. 


Regime III: large amplitude theory (for p g > 0.008.) - For this region 

there are very few available methods of solution. The analysis should incor- 
porate large rotation effects. Equation (3.1) is not valid for such analysis. 
For static cases Reissner's theory has been used with some regularity. The 
problem is a fully three-dimensional elasticity problem with all components 
of motion fully coupled with each other. Little progress has been made in this 
area over the years. Experimental studies in these amplitude ranges are few, 
and most have been using plates which are fairly thick (compared to the thick- 
ness of 6 ym with which we are dealing) . Even before the small deformation 
implied in the theory for Region n becomes invalidated, it happens more often 
than not that the tensile stresses under large amplitude motion go beyond 
tensile strength limits. The structure becomes locally plastic, and we no 
longer can use the isotropy assumptions implied in the analysis. 

For flimsy materials such as we use, the maximum tensile stresses induced 
in large amplitude static motion are given by 


a 


t = 



( 2lV Eh 

W 8 12(1 -v 2 ) 




w 2 1_ 

L 2 8 


+ 


h w 
Z l 


(3.4) 
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Ttiis stress should be below the tensile strength of the material we use, 
if assumptions of isotropy are to be valid- For the periodic structure we are 
using in our analysis, the longitudinal inertial effects are negligible as will 
be shown later in the analysis, and hence the large amplitude vibration problem can 
be viewed as a quasi-steady problem, especially while evaluating stresses. Hence 
the bound on the limits of an analysis such as that for Region II will definitely 
fail if 


o 

where ts = tensile strength of the structure. 

For the polyester film we use as the structural material the tensile strength 
a ts = 1.5 x lo 7 N/m 2 - and E = 4 x 10 8 N/m 2 ; the smallest wavelength used was 
X = 1.814 mm and the thinnest membrane had a thickness of h = 2.54 ym. From 
equation (3.4) we obtain the maximum limit on w/h beyond which the structure 
will be at least locally plastic as 


w 

h 



Ts*4* (1 - V 2 ) 
E x tt 2 *y 2 


_2_ / (1 -v 2 ) g Ts 

TTy y E 


84.0 


(3.5) 


We had given the limit for validity of analysis for Region II by consideration 
of the small deformation approximation implied in it as 


w 

h 


Q. 03 
Y 


21.5 


(3.6) 


for a = 1.4 x 10 3 , a typical case. Hence, it is feasible using the structure 
we are designing to study the limits of validity of the moderate large amplitude 
analysis for the first time in our knowledge. 
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3.4. Structural Response Theory for (Intermediate) 

Large Responses of the Periodic Structure 

Equation (3.3) in nondimens ional form represents the governing equation for 
transverse motion of the structure. We make assumptions (i) , (ii) 5 and (iii) 
of section 3.1 in our analysis; i.e. b >> & and 


ff w dx dy = 0 
v 

furthermore , we assume that the excitation field is of the form 


f p sin dx = 0 for m = odd numbers (3.7) 

J 0 e 1 

where p £ is the loading function in each bay. The excitation field that will 
be derived in the next section can be shown to be of the above form for a given 
form of voltages at terminals A and B of the electrostatic setup. 

For a simply supported rectangular bay, the deflection shape can therefore 
be only of the form 


w 


-22: 


w 


m n mn 


sin 


2mirx 

i 


sin. 


BIZ 

b 


(3.8) 


Furthermore, because of uniformity of 
the choice of the deflection shape is 

w = 221 w sin sin (2n - 

m n mn l 


the electrostatic field in the y 
restricted to 



direction , 


(3.9) 


or, in nondimensional form. 



m n mn 


sin(2irnrx) sin(2n - 1) (fry 8) 


where 8 = 1 / b (3.10) 

defining : 

m = 2mir (3.11) 

(cont 'd) 
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yields 


g- = (2n - 1) gir 
n 


a =ZE 


(t ) . - - . - 

a sm m x sin g-y 
mn n J 


(3.11) 

(concl'd) 


3.4.1. Midplane forces. - During deflection under transverse load the 

midplane suffers in plane strains E , E E Assuming that the external 

xx yy , xy 

load parameter extends through Region II, we define the midplane forces as 


N 

x 


N 

y 


Eh 

(1 - V 2 ) 
Eh 

(1 - V 2 ) 



+ 


+ 


v E ] 

yy/ 

v E ) 

XX / 


N 

xy 


Eh 

2(1 + v) 



) 


(3.12) 


In terms of the nondimensionalization adopted in equation (3.2): 
N = E — + v E— 

x i - v 2 L “ 

Y ± . v 2 l YY xxj 

N = — E — 

xy 2 (1 + v) xy 


(3.13) 


The midplane strains are related to the components of deformation as 


E— 

xx 



1 

2 




8v _1 / 9a\ 2 

3y ^ \3y/ 


E— 

xy 


3u + 3v 2s. IS. 

3y 3x 3x 3y 


(3.14) 
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where u = u/&, v = v/£ and u, v are two dimensional in-plane displacements. 
Before setting up the dynamic equations of in-plane motion we introduce the 
following consideration: 

The frequency of excitation for the electrostatic wall is well below the 
longitudinal natural frequency of the structure. For this case the in-plane 
motions are stiffness dominated. Hence, the effect of longitudinal inertia can 
be neglected. 

The above requirement is satisfied if the reduced frequency fL/c « 1 

fL 

(but a more relaxed criterion can be - — < 0.10) , where f is the operating 

frequency, c is the shear wave speect, and L is the wavelength of the 

primary motion. For the electrostatic model L = i, the maximum frequency of 

operation f = 10 kHz, and c^ = ^90 m/sec. For the largest wavelength 

model that we considered, 1 = 7 mm. Hence, = 7° = 0 04 < 0.1 Thus the 

C sh 1690 

neglect of longitudinal inertia in all the subsequent analyses is well justified. 
The dynamic equations of in-plane motion are now given as 



3x 3y 


3N- 3N — 

—L + = 

3y 3x 


0 


(3.15) 


Substituting the relations given in equations (3.13) and (3.14) we get 


3 2 u , 1 - v 3 2 u . 1 + v 3 2 v . _ ,- - r, 

+ — x + — - + F- (x,y,t) = 0 

3x 2 2 3y 2 2 3x3y X 


3 2 v , 1 - v 3 2 v 1 + v 3 2 u . _ , r, 

— + — 17 + — 2 ITT + F-(x,y,t) = 0 

3y 2 3x 2 3x3y ^ 


(3.16) 


where 
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3s r 9 2 a X - v 3 2 al 1 + v 3a 3 2 a 
F-(x,y,t) = — — — + — — + — 5 7 -7—7 

3x L3x 2 3y 2 J 3y 9x3y 

3a T 3 2 a . 1 - v 3 2 a1 . 1 + v 3a 3 2 a 
F-(x,y,t) =— — + — r 7- + — 7-7^; 

y 3y l_3y 2 3x 2 J 3x 3x3y (3.17) 

For the deflection shape, we use equation (3.11). Substituting in equation (3.17) 
we get expressions for F- and F- as: 

F x = iEEZZ a ik a jA { sin L (I + cos [ (3 k + e i )y ] 

• D !i 2 + 6 k + 3 k 3 i}] + sin [ (I " 5) *] 

• cos[(3 s + 0*)y][j {i 2 + + H 21 6 k ^ 

- Sin [(I + j)x] cos[( 3£ - 0s;)y][j{l 2 + ej - B k 3 i)] 

- sin [(I - j)x] cos [( 6 - - 0 j)y][j{? + 2 k 3 ^)] (3.18) 

F y - jEEZE a ik a j* { [( 0 * + 0 j)?] cos [(I + j)x] 

• |]0j(0| + i 2 + 1 2 ^ i 3 )J “ sin £(0£ + 0^)yj cos [(i - j)x] 

• [ 6 J <6 I + 12 - HP 1 5) ] * sin [< 6 fc - 65 )?] [li + 5)S 

' [ S Z (6 i * ~ i V + 1 2 V i 5>] - sin [(Bj - 8 j)y] cos [(I - j)x] 

• (3.19, 
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From the known form of F- and F- one is now able to proceed to 
solve u, v. 

3.4.2. Edge conditions. - The edge conditions describe the state of 
fixity of the structure for in-plane motions along the edges. For the simply 
supported periodic structure the edge restraint is 

u(o,y,x) = u(l,y,T) = 0 

v(x,o,x) = v(x,j,t) = 0 (3.20) 

3.5. Solution of the In-Plane Motions 

Since equations (3.16) are time-dependent equations, a transient solution 
which depends on the initial condition will exist along with the solution of 
equations (3.16). This particular solution will be representative of the initial 
displacement fields that may exist (e.g., if initially the membrane is stretched 
with uniform tension T, the initial displacement field will be nontrivial) . 

For the case of zero initial tension, the particular solution takes the 
trivial form: 

u o (x,y,T) = v Q (x,y,T) = 0 (3.21) 

For the solution of equation 3.16 we take the displacement fields as. 


rzrz [\ jia 

(t ) sin | (i + j ) x | cos 

K + 


+ B ijkJi <T) sin 

| (i - j ) x | cos | (6j + 



+ c ijM (T) sin 

| (i + j ) x | cos ( (6- - 



+ D ij k » (T) sin 

{ (i - j ) x | cos { (6- - 

^0 

(3.22) 
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-ZE£Z[E ijkJ <t) c°sja+ j)j)sin{(8j +8j>y) 

+ F ijk» (T) cos ( (I ■ 5>k} s p< B s +B i>?) 


+ s ijk* (T> C H (I + ?>*■} s H (B k - 8 j>? } 


H ijki (T) “*! (I ■ 5,i l sin ! (B k 


(3.23) 


Substituting for the derivatives of u, v, F-, F- in equation (3.16) and 
equating like terms one obtains a set of algebraic equations for the coefficients 
A ijkJl' etc * We 9^- ve one set these equations as: 

A ijk* ( l + i) + + 6 i) J + E ijkji [Hr^ ( l + ^)( 6 k + 6 a). 


a. 3^. 

ijkJl 2 


= I a ik a j*[ J | 12 + B | + P^ B k B l)] 

P ( B k + B k)( I + j)] + E i jM [( B k + 6 l) * P^ ( J + J ) ] 

= 4 a ik a j)t[ B i I B ! + P^ * 2 + P^ 1 i). 


(3.24) 


(3.25) 


and similar equations connecting other coefficients. Finally the coefficients 
are given as: 


A ijk& 


a ik a j £ 1 

2(1 "’{(I + j) 2 + <6- - ftj) 2 


_5 { iP - 


P B S + HP B k P ( B k + B i) 


1 - V /- 


(; * j) 


•(PP ( J + J )(«k + P b j(4 + PP * 2 + PP 1 0. 


(3.26) 
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For i J- j 


B . = + 


a a . „ 
lk ]£ 


ijk& 2(1 - v) 


L(i - D ) 2 + ( e - + Zj ) 2 } 


• ^ + 4 + b e ej {(b e + B-,) 2 + ip(l- 5) : 

H)] 


( l - 5)( 6 e + ' 6 i) e i (4 + - 


- v -*■ 1 + V 

2 1 " 2 


(3.27) 


when i = j, B . , . „ = 0 . 
J ' :iok& 


C. 


a., a. . 
lk 3& 


ijk£ 2(1 - v) 


L(I + j ) 2 + (3- - 3j ) 2 


4 • B k 6 i][( s k - 6 j ) 2 + H-K 1 + j) 2 ] 


pu)} 


+ ip (i + j)(sj - Bj) Bj (e| 


1 - V T 2 . 1 


1 + 


(3.28) 


and whenever i - j, D. = G, and in other cases 
J 13 k£ 


D. 


a., a. 
3- k 3*- 


ijkji 2 (1 - v) 


-.2 


L Cl - j ) 2 + - 3 j) 2 J 


• [ j (i + a k a l) ( a k + a I + H" 2 1 ~ ^ ) 


.^(i . j) B . ( Bj . - 6 j)(b| 


1 - Vt 2 1 + Vtt 


1 - 


u)] 


(3.29) 
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2 


a., a 
ik 

2L 

1 




E ijk£ 2(1 - 

v) 

_(i + j) 2 + 

+ 8 i )2 -i 



•[^ 

■fe 

+ e i) ( L + 3 

Ml 1 ' + i 

• r 4 - 

f 2 6 k 

•{(M 

2 

1 + 


8 s) ) ! 4 

+ H- 

1 + 2 1 ^ 1 e i] (3.30) 

a., a 

p lk 

2L- 

1 

- 

2 


ijkA 2(1 - 

v) 

_(I - j) 2 + 

(8j + 6j> 2 . 





+ e i) ( l - $ 

j IV + ±- 

V4 * 

1 + v 1 

2 ^k 

•[»■ 

5 ) 2 

+ 2 ^k 

+ s i) 2 ] [ 6 1 


(3.31) 


3.,,. = 0 for k = A and 
13 k£ 


ijkf, 2(1 


a ik a j& 1 

2(1 " V) L(I + j) 2 + (gr - g-) 2 _ 


k w l' 


{H^ k - »i) (* + 5) 5[i 2 ♦ H* 4 - H* »s »r] 

■ <•$ ^T^k- 8 *-) 2 ] [4 + 1 ?] ) <3. 


32) 


and H..,„ = 0 if k = £ or else 
xjkil 


H. . 


a a 
ikji 


ijkA 2 (1 - v) 


L(i - j) 2 


- i il 

+ (6 k - ^ )2 -l 


'(H^ ( B k ■ e l) ‘ 5) j [i 


A - I - 2 1 - v 2 

+ — 8 ; 


1 + v 


3 k j 


- ** |> - ^) 2 + k - fi-tft + ^i 2 - 3]} 


i j I > (3.33) 
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Having obtained the coefficients one is now able to evaluate the midplane 
forces through the help of equations (3.13) and (3.14). Thus 


N = 
x 


ZXZX k + j ) + V E ij kJi (^ k + 


(1 - v*) 


a i3c a j£ A t 


(l j + v 3- 3j) } cos(l + j)x cos ^ + Bj): 


’ ( B ijk& ( x “ + v F ijkJl ( 0 k + 6 l) 


a., a . ( 

ik -ii 


i j - v S; 


S >l)] 


• cos 


(l - j)i cos(8 s + Sj)? 

kt l 1 + + v G ijM ( 6 ic ” 6 i) 


a ik a j£ / T t 


(l j - v S E Bj)} 


(i + j)x cos (bj - Bj)j + j =03 (i - j); cos (e- - ej)y) 


+ “ijki 


(1-5) 


+ V H. ... 
ljkJl 


k ■ 8 -J 


ik jZ (r t 
+ g- 1 - : + V 3 


J 6 i)| 


^—ZEEE [{v (i + 5) + E i 3M (»s + 6 s) 


(1 - \» ) 


a., a 
ik ]£ 


(sj 6J + V i j) } oos(l + j); =os(b e + B^i 


f B ijkB M) + F ijk* ( 6 k + ^ + ^ (s k B t - i 5 v)) 
oos(i - j); oos(b s + Bj)? + {v c ijkk (l + j) + G. jkk (sj - B^ 


a., a. „ 
ik ]£ 


[v i j - B k Sx]Jcoa(i + 5)5 cos(s k - B^i 


(3.35) 
(cont ’d) 
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♦ (V D_. (i - j) + H. jM (6jj - Sj) + [v i I +% % 


"ijkA 


cos 


and 


N 


xy 2 (1 


(i - j)j COS (l5j - B ? )jj 

{% j) - A ijkt (ejj - Bj) 


- E «W ( Z + j) } sin ( J + *)* sin (?k + 8 l) ? 

(»k +8 i) + F ijkn - J ) + -^T 21 8 k 5 } sin ' j) 
sln ( 8 k + + (" C ijkji ( e k " 8 a) " G ijka l 1 + 


+ 8 k 5 } sin ( I + i)* sin («k - 8 l) ? - { D ijk 4 ( 8 E - 8 l) 

( z - 5 ) + B k J - 5)* Sin ( s k " s ^) 5 ] 


+ H. M 
ljkJl 


Having obtained the terms N^, N^, one can now solve equation (3.1) 


3.6. Solution of the Transverse Motion 

The governing equation (3.3) can be solved as modal equations in time. 
The form of a is given by equation (3.11) as: 

* -EE a sin m x sin(3n y) 
m n 

The Raleigh-Ritz approximation to equation (3.3) is obtained by multiplying 
equation (3.3) by sin p x sin 3“ Y and integrating over the x - y space. 


(3.35) 
(concl ’d) 


(3.36) 
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Defining: 


and 


and 


Jl) 

ijk£ 


A (2 > 

13 


a (3) 

ijU 


A (4) 

ijU 


B (1) 

ijk£ 


b {2) 

B ijkJl 


b (3) 

ijki 


b (4) 

B ijk& 


c. {1) 

i jkS, 


c. 


( 2 ) 


ijk 


[M 

A ijkf 

+ v E ijU ( 

?k + ^ 

£ B ijkf 

M 

) + v F ijkt 

(«k + 8 j) 

[ C ijkJl 



( 6 k ' 8 i) 

[ D ij k 2. 


) * v H ijk 4 

(?k - ^ 

V 

( I+ 

5 ) + ( 

'& * *i) 

V B ijkf 

(f- 

5 ) + F ijki ( 

k + e s) 

v C ijk* 

(- 

i) + ^ijki ( 

fe - 6 i) 

v D ijk& 

(- 

i) + H ijki ( 

fc - e s) 


i Vil e i: 5 ‘ ( 6 k + B i) - 

- [ B i 3k s. k + 6 i) + ( r - *) 


a., a . 

ajc 2 & 
8 


a. . a 
8 


a., a, 

+ 

8 


a., a. . 
lk jS, 


a., a . „ 
ik 

8 


a., a 
lk i* 


a a 
ik M 


a. , a . . 
ik 3 A 


(i J + V B E Sj)] 

(i 5 - v 8j)] 

(i J - v BjJ Bj)J 

(U - v 8 E S^)J 

(Bj; Sj + v I j) 

( 6 s 8 i - 1 i v ) 

( v 1 5 - 8 k e i) 

( v 1 5 + 6 k e «) 



(3.37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 
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(3.47) 


C ijk£ -C ijkJl ( P k " & j) 

and 

C ijk£ = _ [ D ijk£ ( e k “ B l) 
and also the following parameters 
Hi - i + j + p 
R 2 = i + j + r 
R 3 = i + j - p 
R 4 = i + p - j 
R 5 = i + r - j 
Rg = i + p + r 
R 7 = j + p + r 


G ijk» ( l + 5 ) + 

+ H i jM $ - i) 


a M a * 0 

B- j 
4 p k J 


a., a . . 
lk ]£ „ 

* 4 


'] 


and 


51 = k + Jl + q 

52 - k + i, + s 

53 = k + s + q 

54 = & + s + q 

55 = k + l - q 


(3.48) 


(3.49) 


(3.50) 

(cont'd) 
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S 6 * k - i + s 


S 7 - k - s + q 

and the Kronecker delta function 
£ 

rs = 0 r ? s 

5 

rs = 1 r = 1 


the Rayleigh-Ritz approximation equation (3.3) is then 


d 2 a 


JB3. 


+ K a 


dt 


2 P<2 P<1 


+ mzumm — - — a a ik a i 0 

i j k A r s 4(1 - v 2 ) pq ik 


[ 


$1 Pi + 62 P2 + <$3 ^3 + ^4 P4I- P 




pq 

'Gen 


where 


K = ll /-? ■ ~2\ 2 

Pq 12(1 - V 2 ) 


♦ $ 


p i = — - — fp 2 a.^?_ + e- e. (1 ?. 

- 1 a, .a, „ ljki q i]ki 


ij kJL 


p 2 = 


P 3 = 


P4 = 


iiV t 


Z 2 . ,(2) q 2 q (2) 

a ij a ki |P ijk5, 8ijk£ 


2) 1 

)k£j 


^ [? *SL • i «.] 

ij si [ pJ ‘ e 4*] 


and 


(3.50) 
(concl 'd) 


(3.51) 


(3.52) 


(3.53) 


(3.54) 
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6i = U + 6 - ) (d + 5„ - 6 o ) 

1 \ Rlr R 2 p R 3 r / \ Sis S 2 q S 5 S/ 

<5 9 = ($ +6 - 6 ) (6 +5 -5.-6.) 

02 \ Sis S 2 q S 5 s J \ R 4 r R 5 P Rej R7V 

63 = (5 +6 -6 ^ /6 + 6 - - 6 o N 

3 V Rir R2P R3 r / \ S 7 £ S 6 q S 3 & S 4 k/ 

64 m U + 6^ - s n . - 5 \ ( s e + S B ~ 6 - 6 \ (3.55) 

4 \ R 4 r Rsp R 63 R 71 / V S7£ S6q S 3 JI S4k/ 

and 

P Gen = ^ X P e sin x ^ sin ^ dx dy J (3.56) 

The set of ordinary differential equations is highly nonlinear because of the 

3. cl 3. ' 

presence of terms such as ij k£ pq. The coupled set of such equations can 
be solved using numerical techniques. Analytical solutions might be impossible 
except in the simplest case, e.g. a one mode solution such asi=j=k=£= 
p = q = 1. When the electrostatic wall model is subjected to fluidic forces such 
as will occur if used for drag reduction studies, a complete analysis of the 
problem as formulated above with a few normal modes in the deflection shape might 
be necessary. This is especially so in flutter regimes of structural motion. In 
such cases the generalized loading will consist of three main contributions : 

(a) random turbulent pressure loading which can be insignificant unless flow speeds 
are fairly large and the membrane thickness very small; (b) wall pressure loading 
due to the interaction of a pulsating boundary in a turbulent boundary layer flow; 
the magnitude of this loading can be significant under a variety of circumstances 
including the flutter mode instability case; and (c) the electrostatic wall loading 
due to the primary excitation through electrostatic attraction forces. 

In an actual experiment with flow over the structure such as in a windtunnel 
other loadings might also be present. The most common of these are pressure 
gradients in tne tunnel and static pressure differentials across the membrane. An 
accurate prediction of the structural response depends largely on identifying all 
these influences. 
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4 . ANALYSIS OF THE ELECTRIC FIELD 


4.1. Determination of the Primary Electric Field 

The schematic of the electrostatic system is shown in figure la. The output 
from the transformer T is biased at the center tap and connected to terminals A 
and B as indicated. An array of electrodes which is etched on a printed circuit 
board and coated by degassed epoxy cement or lacquer forms one part of the 
electrostatic wall system. A grounded conducting sheet (membrane) which is struc- 
turally supported periodically at distances apart and separated from the 
electrodes at a height "H" forms with these electrodes a capacitive network. 

The membrane is free to move transversely? when alternate electrodes are connected 
to the terminals A and B respectively the membrane is subjected to an electro- 
static force distribution and deflects under these forces. The space between 
the membrane and the electrodes is filled with air and sealed. To obtain maximum 
force for given terminal voltages, the distance between the membrane and the 
electrodes should be as small as possible. However, as the distance between the 
membrane and the electrodes becomes smaller and smaller the applied voltages in 
the terminals must be reduced lest a breakdown of the electrostatic wall system 
due to arcing will occur. The breakdown potential for a given separation in air 
is given by Paschen’s law. For air at NTP the breakdown rms potential gradient 

is 3.1 kV/mm. If instead of air a fluid such as SFg is used, the breakdown voltage 
gradient can be raised by a factor of two. 

In figure 2 we show a typical bay of the electrostatic wall configuration. 

The electric field between the electrodes and the membrane can be obtained by 
solving in the domain 

V 2 V = 0 (4.1) 

with the boundary conditions (as indicated in fig. 2) . 

z = 0 V = 0 

z = H V = Vq + Vi (x,H) sin fit 

x = 0,1 V = Vo z/H (4.2) 

In figure 3, Vi (x,H) and its Fourier series representation is shown. 

From the figure it can be seen that 
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(4.3) 


Vi(x,H) * 1.1463 Vi sin 27r-r- sin Qt 


is a good approximation. 

The solution of equation (4.1) with the boundary conditions given in equation 
(4.2) is rather straightforward and can be written as 


V z 

o 2 ttv 

V = — — + 1.1463 V ]_ sin — r— sin Qt 

n 


sinh 




sinh 2tt ^ 

Z J 


(4.4) 


The electric field in the dielectric between the electrodes and the membrane 
is thus 


E = E - O 
X y 

V Vi cosh [4^3 

E ~ ~~ + 7.2024 — sin sin Qt — 

Z H jL Jv 

sinh 



(4.5) 


Force density due to primary field. 



P o 2 z — 0 


The membrane is subjected to a force 

(4.6) 


Thus the generalized electrostatic force is 


p£ q =0 for p > 1 
Gen * 


and for P = 1 


pPq . 

Cen (2q - 1) 


sinh 


(¥) 


'0. 


935 e v v, 
o 1 

EyHJl 


■] 


sm to t 


(4.7) 


Optimal values for V q , V, . There is a limit to the maximum voltages 
Vq, V permissible for a given configuration, which is determined by the rms 
voltage gradient to be less than the breakdown value (3.1 kV/mm for air at NTP) . 
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From equation (4.5) the maximum rms electric field value is given as 



<_ 3.1 Mv/m (4.8) 

In order to obtain the optimal values of Vg, Vj for maximum generalized 
forces, one must minimize the generalized forces subject to the constraint 
equation (4.8) . 


4.2. Requirements, Ratings, etc. of the Electrical 
Networks for Electrostatic Wall Configuration 

In figure 4a is shown the schematic of the electrical setup. Power is 
drawn from a signal generator stepped through a power amplifier and fed into a 
transformer. The transformer output is fed into the exciter for the electrostatic 
configuration. Basically the electrostatic wall comprising the membrane and the 
terminals act as a capacitive load. The capacitance can be measured using a 
capacity meter available commercially. It is also fairly straightforward to 
theoretically model this electrostatic wall as a discrete capacitive network and 
evaluate the capacitances of the wall as done in reference 6. 

For satisfactory performance of the electrostatic wall system the following 
points should be kept in mind: 

(i) The electrostatic wall should be operated below the breakdown voltage 
levels to avoid arcing and sparking and consequent degradation of the terminals 
and burnout or charring of the membrane surface. 

(ii) Collapse of the membrane into the cavity and consequent shorting of 
the electrical system should at all costs be avoided lest the transformer or 
power amplifier be damaged. 

(iii) Any L-C oscillation due to the loading of the system should be avoided. 
This is accomplished by designing the outages of the units such that the external 
load is well within the operating load for the system. 
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4.3. Evaluation of the Maximum Permissible Capacitive 
Load for the Electrical System 

The ratings of the transformer are 

Turns ratio a = N 2 /N 1 = 40 
Resistance on the primary side Ri = 0.5 J2 
Inductive reactance on the primary side Xi = 8 ft @ 1 kHz 
. Resistance on the secondary side R 2 = 2 kfi 
Inductive reactance on the secondary side X 2 = 50 kQ @1 kHz 

External load is capacitive with a reactance X c @1 kHz = t~ 

In figure 4b the equivalent circuit is shown. It is assumed that the 

magnetizing current is zero, an assumption which is quite valid (usually it is 

5 to 10 percent of the primary current) . With this assumption the reactance 

X = co. 
mag 

From figure 4: 

The equivalent impedance 
primary side = Zeq 


The primary current 
The load current 

The circuit becomes a resonant circuit when the reactance is zero, 

i.e. X T + X T = — — (4.12) 

L 1 J* u qa 2 

a 2 

The operating range of the electrostatic wall for experiments lies between 
200 Hz to 2 kHz« In order to protect the power amplifier the resonant frequency 


on the 


l 2 r 2 1 

= j X r + Ri + j + — + 

1 a 2 a 2 jojCa 2 

X. 


- tx + — 
L 1 a 2 


wCa'* 


j - R! + 


a‘ 


= 


II = 


V 

eq z 


v/z_ = i L -e 


(4.9) 

(4.10) 


I 2 = ll/a = — L -Q 
za 


(4.11) 
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which can be determined from equation (4.12) should be well above this range. 
Fixing the resonant frequency f to be around 3000 Hz we obtain the condition 
that capacitance of the electrostatic wall be 


c < 


= 280 pf 


_ . Xlf 
27rf 1000 


a 2 + 


X2f 

1000 


(4.13) 


With a proper choice of equipment (i.e. low reactance transformers and 
power amplifiers) , the size and capacity of the electrostatic wall can be 
enhanced considerably for the same range of operating frequencies. The fore- 
going discussion makes it imperative that when different models require 
laboratory testing the capacitance of the model be measured to see whether 
there is any limitation on operating frequency imposed by the choice and 
availability of power equipments. 

5. THEORETICAL ANALYSIS AND EXPERIMENTAL VERIFICATION OF THE 
PERFORMANCE OF THE ELECTROSTATIC WALL 

5.1. Analysis of Structural Motion 

We consider the case where the generalized force is due primarily. to the 

electrostatic field set up in the model. The structural motion resulting can 

be studied by solving equation (3.52) with the expression for provided 

Gen 

by equation (4.7). The deflection shape a is given in equation (3.11) as 

a a sin m x sin 0- y 

m n mn n 


We consider the case where 0 = £/b << 1, corresponding to the situation of the 
width of the bays b - 50 l which was true for most of the wall to be built for 
tunnel testing. 

From the nature of the generalized load P^[ given in equation (4.7), 

i.e. 


P 


pq = 

Gen 


0 p > 1 
F (q) p = 1 
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it is obvious that the deflection shape can be approximated as 


a 



sin 2ttx sin g- y 
n 


(5.1) 


It is not at once clear how far the summation Implied in equation (5,1) must 

be carried out. However, the nature of the generalized forces P^ (being 

1 Gen 

proportional to ^ suggests that at best, the first two modes are important. 

The equation system for two mode series is of the form 


d 2 an 

— + ail 

dx 2 


M 


+ G 


11 a ll 


+ G 


12 a ll a 12 


+ G 


13 a ll a 12 


P 


11 

Gen 


d 2 a 

f- + a 22 [k 2 2] + G\ z a i2 + G 2 1 a ll a 12 + G 23 a ll a 12 = P qL (5 ' 2) 

In table 1 the coefficients g}^ G^f g 13^ kill etc. are all tabulated 
for arbitrary 3/ using a program called MACSYMA available through ARPA, NBS 
or MIT network. 

5.1.1. Initial approximations . - We shall examine the case of 6 << 1 
and assume that the percentage of the second normal mode in the solution to the 
problem is entirely insignificant in order to obtain an idea of the sensitivity 
of the solution to truncation of the series. Setting 3=0 and taking only the 
first of the equation set in equation (5.2) one obtains 

d 2 au 

— + *11 ail + S 1 !! a n - p“ n (5.3) 

dx^ 


From equation 
where Pig = 


(4.7) 

1 

sinh 2irH 

l 


p}, 1 = Pin sin fix 

Gen iU 

0.935 e V 0 Vi 

iyil 


(5.4a) 

(5.4b) 


Then terms kn and in equation (5.3) have (for 3=0) values 
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k ll = 


16tt 4 y 2 


(5.4c) 


12(1 - v 2 ) 



(3 - v 2 )tt 4 
(1 - v 2 ) 


Equation (5.3) can be rewritten using the following transformation: 



and ?2 = E*io 



(5.5) 


as 


d 2 u 

d0 2 


+ u + u^ = P 2 cos 9 


(5.6) 


5.1 # 2. Harmonic solution of equation (5.6). - Seeking harmonic solutions 
for u, equation (5.6) leads to a simplified cubic equation. 


(3 i - 1) u + ~ u 3 = P 2 


(5.7) 


from which u, and hence a, can be obtained. 

5.2. Theoretical and Experimental Results 

Theoretical results for two configurations of the electrostatic wall are 
presented in figures 5 through 8. The membrane material used (commercially 
named Mylar) has the following properties: 

E = 3.5 x 10 8 kgf/m 2 

p = 138 *S£.3 g .£ 2 

m 4 

v = 0- 3 (5.8) 
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The air gap between the membranes in two different configurations were 
H = 0,305 mm and H = 0,127 mm. Since the electrostatic forces are inversely 
proportional to gap height H, the smaller gap provides a larger force on 
the structure for a given voltage. However, the breakdown potential for the 
system is smaller when the gap is smaller. Hence the electrostatic forces 
available in both these configurations at the optimal voltage are of the same 
order. However the smaller air gap has the advantage that one can work with 
lower voltage levels. 

In figures 5 and 6 the theoretical values are compared to the bench test 
results. It is surprising to observe such an excellent agreement between the 
test values and theoretical prediction. The experimental measurements of 
dynamic surface motion were obtained using an optical setup. Reference 6 
describes in detail the measurement technique. The a.c. voltage applied to 
the electrodes was at 300 Hz for these experiments. 

Figures 7 and 8 show the frequency response of the structure at various 
excitation levels. The nature of the backbone curve is indicated by dotted 
lines in these figures . At low excitation levels where w/h « 1 , there is a 
steep increase in amplitude levels at resonance. At higher levels of excitation 
resonance has little effect on amplitude level. In the theoretical analysis 
structural damping was taken to be negligible. Since damping primarily affects 
the near resonance amplitude levels, it is not necessary to include structural 
damping in the response studies for the present case as can be seen from the 
nature of the response curves. The reason for this is very clear? viz, the 
test structures have resonance frequencies well above the excitation frequency 
for load values of interest (large amplitude case) . Furthermore at these load 
values the nonlinear stiffnesses themselves act as a damper or delimiter on 
the amplitude levels. When the structural nonlinearity is of the soft spring type 
it might well be important to include the damping in the analysis since natural 
frequencies for this case are being constantly shifted downward from the low 
excitation case (w/h << 1) with larger excitation levels. 

5.3. Examination of the Perturbation Field 
Due to Primary Motion 

When the membrane is set into motion by the action of an electrostatic 
force field, the electrical field configuration changes and consequently the 
excitation field itself must be reexamined. 
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To analyze the effect of the surface motion on the electric field configu- 
ration an orthogonal curvilinear coordinate system can be used- Since, as 
far as the electric field is concerned, the structural motions have a station- 
ary behavior, a quasi-rigid coordinate mapping is required- Curvilinear 
coordinates such as the one to be described have been used by Benjamin (ref. 7) 
and Chang (ref. 8) for flow problems- The first order curvilinear orthogonal 
map is shown in figure 9 . The amplitude of the wavy surface is taken to 

be (0.1 H) ; Yj = 0.0 curve in the figure represents the approximation to the 

£ 

actual surface. Since 77 « 1, the variation in the chordwise directions have 

D 

been neglected just as in section 4. The mappings for this configuration are 
xi = x - a e" 2lTZ cos 2 ttx sin Qt 

yi = 2 - a e“ 2Trz sin 2irx sin fit (5.9) 

x - y coordinates can be written as 

x = xi + a e -2Tr ^l cos 27rxi sin fit 


z - Vl + a e sin 2 ttxi sin fit 


(5.10) 


The solution of equation (4.1) with the boundary conditions should be attempted. 
We write the governing equation and the boundary conditions in the mapped system 
as 


3 2 V a 2 v 

+ = 0 

3xi 2 3yi 2 


(5.11) 


Yl = 0 V = 0 

yi = (H/£) v = V Q + vf sin fit sin(2irx 1 ) 

xi = 0, 1 V = V qY i (5.12) 
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The solution of the above problem yields 


V 


z 




sinh 


2ttH 

l 


sin 2 ttx sin fit 


+ 27ra(cos 


2.; 


- cos 4mc) sin fix 


(5.13) 


and 


V = 
x 


0 • 


The correction to the electric field distribution is therefore of the order of 
(w/£) times its original value. For the electrostatic walls we are testing 
the magnitude of this correction is less than five percent when the amplitude 
of the surface motion is 0.025 mm at a spacing of the bay of 1.814 mm. In all 

the tests that we carried out the perturbation field was negligible since the 
excitation field could barely drive the system to these values of amplitude. 

An error of 15 percent in the evaluation of the electric field would usually 
generate an error of 10 percent in the amplitude prediction, which must be 
borne in mind when accuracies of that order are required. 

When the deflections of the surface are quite large the perturbed field 
can be evaluated using a more refined (second order) coordinate system such as 
used by Chang (ref. 8) . We show the mapped system in figure 10 and give the 

c..'V 

transformations below as 

xi = x - a e~ 2lTZ cos 2irx sin fit - ira 2 e _1+7rz sin 4irx sin fit 

yi = z - a e 27Tz sin 2irx sin fit + TTa 2 (e -lf1TZ cos 4ttx - 1^ (5.14) 

x = xi + a e -2 ^ 1 cos 2irxi - ira 2 e _lt1T Yl s i n 4 ttx7 

z = yi + a e 2iry l sin 2irxi + ua 2 ^e -4lT Y 1 cos 4irxi + l) (5.15) 
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5.4. Inclusion of Aerodynamic Forces for 
the Electrostatic Wall 

The electrostatic walls were purposefully designed as small wavelength 
configurations. Hence, static divergence will not be achieved in any of the 
electrostatic walls under operating condition (15 m/s - 50 m/s) in the 7 in. x 
11 in. tunnel. In reference 1, an in viscid analysis has been carried out for 
the evaluation of the aerodynamical ly induced pressure p [eq. (39) , p. 16 of 
ref. 1] due to structural motion. 

We give this expression as 

U 2 

aero » - - 

p = 2irp<air> ai i sin 2irx sin ny sin wx. (5.16) 

Vl - M 2 

' CO 


Figure 11 shows the effect of including aerodynamic loads in the analysis 
for the case of flow over the membrane. Our contention is fairly clear; there is 
no need for any sophisticated analysis of aerodynamic effects for these configu- 
rations since even the inclusion of the inviscid values which overestimate 
these magnitudes increases the levels of amplitude very slightly and that increase 
is only over a very narrow window. The a.c. voltage applied to the electrodes 
in this example is at 300 Hz. If the' a.c. voltage was applied at very much higher 
frequencies the effect of the aerodynamic load on the motion will be enhanced 
since the effective nonlinear stiffness would be considerably lowered. 

We again point out that the analysis which took into consideration the 
aerodynamic load was carried out under the prior known fact that the wavelength 
of the structures was much too short to cause static divergence; hence the flow- 
structure interaction problem is not of an eigenvalue type but merely a forced 
response problem. 


6 . CONCLUDING REMARKS 

A unified theory for an electrostatically driven active wall system has been 
presented. The electrostatic wall system is capable of producing deflections of 

many orders of wall thicknesses. Consequently a large intermediate response theory 
has been used for analysis. The theoretical analyses are compared with bench test 
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results which show excellent agreements between the two. The case of an elec- 
trostatic wall vibrating in the windtunnel is considered, and the theoretical 
predictions under simulated flow conditions indicate that the aerodynamic 
effects are negligible. The perturbation effects to the electric field due to 
the structural motions are also considered and shown to be negligible for the 
test experiments* 
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Table 1. Stiffness coefficients. 
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Figure 2. Mathematical representation of the electrostatic field problem* 



SINE APPROXIMATION 
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Figure 3. Vi (x,H) and its Fourier representation. 


Equivalent circuit 


Figure 4. Schematic of the hookup for the electrostatic wall system. 
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CASE 1: MEMBRANE THICKNESS = 2.54 pm 



Figure 5. Theoretical response curve and test results. 











CASE 2: MEMBRANE THICKNESS = 6.35 Jim 



Figure 8. Frequency response curve. 
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Figure 9. First order orthogonal coordinate system for analysis of the perturbed field. 
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Figure 10. Second order orthogonal coordinate system for analysis of the perturbed field. 
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